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OPTIMIZATION  OF  AUTO-PILOT  EQUATIONS 
FOR  RAPID  ESTIMATION  OF  HELICOPTER  CONTROL  SETTINGS 


ABSTRACT 

An  automatic  feedback  system,  based  on  continuous  monitoring 
of  control  loads,  is  used  to  find  the  control  settings  that  are 
required  to  obtain  a  given  flight  condition  of  a  helicopter  rotor. 
A  program  is  developed  that  searches  automatically  for  the  optimum 
gains  and  time  constants  of  the  system.  Satisfactory  results  are 
achieved  for  given  conditions  as  an  example. 
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OPTIMIZATION  OF  AUTO-PILOT  EQUATIONS 


FOR  RAPID  ESTIMATION  OF  HELICOPTER  CONTROL  SETTINGS 


1 ._ _ INTRODUCTION 

Previous  work  in  the  area  of  helicopter  stability  and  vibrations 
has  shown  that  an  accurate  knowledge  of  the  helicopter  control  settings 
is  a  necessary  prerequisite  to  the  determination  of  blade  damping  or 
rotor  loads.  The  mathematical  formulation  of  this  problem  involves 
solution  of  a  set  of  non-linear  differential  equations  for  the  periodic, 
equilibrium  solution.  This, in  itself,  is  fairly  straightforward;  but 
the  problem  is  complicated  by  the  fact  that  the  unknown  control  settings 
appear  as  forcing  functions  (  and  sometimes  as  coefficients  )  in  the 
equations.  These  control  settings  must  be  chosen  so  as  to  satisfy 
certain  integral  constraints  on  the  solution,  namely  that  the  helicopter 
be  flying  in  trim  at  the  desired  flight  condition. 

In  reference  1,  a  solution  to  the  above  problem  is  formulated 
whereby  a  set  of  control  equations  (  called  an  "automatic  pilot"  ) 
is  used  to  bring  the  controls  to  appropriate  values  simultaneously 
with  the  solution  of  the  blade  equations.  The  coefficients  of  this 
controller  (two  gains  and  two  time  constants)  are  chosen  by  trial 
and  error  to  give  the  most  rapid  convergence  to  the  desired  settings. 

Although  the  results  in  reference  (l)  are  very  satisfactory, 
there  are  several  aspects  of  the  problem  which  merit  further  study. 

First,  the  choice  of  parameters  in  (l)  was  made  on  the  basis  only  of 
a  qualitative  assessment  of  when  the  controls  had  converged.  A  more 
quantitative  (and  automated)  approach  is  needed  before  the  method 
can  be  extended  to  general  problems.  Second,  several  errors  have 
been  found  in  the  equations  of  motion  of  (1).  Thus,  the  results 
need  to  be  verified  for  an  accurate  set  of  equations.  Third,  the 
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results  in  (l)  do  not  include  a  study  of  convergence  or  optimality 
conditions  (local  minima,  etc.) -Therefore,  a  more  complete  analysis 
is  warranted. 

In  this  paper,  the  corrected  equations  are  studied  in  more 
detail  in  terms  of  convergence  properties.  Whereas  reference  (l) 
concentrates  on  loss  of  stability  in  extreme  conditions  (stall, high 
advance  ratio,  low  torsional  frequency) ,  the  present  work  concentrates 
on  the  controller  characteristics  in  the  normal  operating  range. 

To  this  end,  a  program  is  developed  that  searches  automatically 
for  the  optimum  gains  and  time  constants  of  the  system. 
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2.  MATHEMATICAL  DESCRIPTION 


2.1.  Rotor  Equations 

The  physical  and  mathematical  models  used  here  are  the  same 
as  in  reference  (l)  with  the  exception  that  certain  algebraic  errors 
in  (l)  have  been  corrected.  The  physical  model,  given  in  Figure  1, 
shows  a  single  section  of  a  slender,  rigid,  inelastic  blade,  which 
is  hinged  in  the  torsional  and  out-of-plane  directions  at  the  center 
of  rotation  with  restraints,  Kc  and  Kg  .  The  blade  is  assumed  to 
flap  with  angle  <3,  and  to  feather  with  angle  .  Fixed  coordinates 
«ue  dof:i  nod  with  the  Z-direction  along  the  rotor  shaft  and  with  the 


X-direction  opposite  to  the  direction  of  flight.  The  blade  rotates 
about  the  shaft  in  the  XY  plane  with  constant  angular  velocity  ft  . 
The  rotor  shaft  angle  Y,  which  is  measured  from  down  wind,  is  in 
radians.  The  blade  position  with  respect  to  the  fixed  coordinate 
system  is  thus  defined  by  the  three  angles  \p,  ft  ,  and  0 ,  which 
uniquely  define  the  blade  position. 

As  the  result  of  the  derivation  and  simplication,  we  have 

%  ♦  f/3  =  Fa  r  (p*--  ,)#r  (1) 

8  +  6  =  Mc  -  (  t£;  -  i)  (6  -  -  £  sin  -2t  )  (2) 

The  aerodynamic  force  F^  and  the  moment  about  the  pitch  hinge 

Mfl  are  obtained  from  piecewise  linear,  quasi-steady  strip  theory, 

PA  -  Jf  (;*(  e  ■  Fp  /  Ur  )  r  £  fac'd  L\-  (  P  rJL/2  -  <]>  ) 

Mff  =  ■  'Tz  f  Url(  P  +  -t-d  *  r  T  f  ^ c  d  Ur  C,n 

Combining  (l),(2),(3)i  and  (4),  we  have 

jj  -T  f(  I  t J*.  Stn  T)/2>  r  l pl  t  (  t  r/U!>,nY)(~jUCcj></  - 

-  $  -=  £  i  l.  (  t  f~) 

C  G  -  <t> )  —  (  i  t  /-v  i.Viyv  xj  +■  ( ,)  ^  (.5) 

o  ''h*"***  +  c4l<?  -0( 

-  (  a'/-  l  )(  PL  r  <Ti  5,V  if  -f  £  ccsf  )  r  fj~  ( 

z  Iff-  (  '  f~:  C  I  f  3  f 

,-177.1  I  t^S,nf)^  -r  Co 


(3) 

(4) 
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2,2.  Auto-pilot  Equations 

In  order  to  simulate  the  trimmed  flight  of  a  helicopter  rotor, 
the  rotor  is  maintained  at  a  fixed  value  of  the  thrust  coefficient 
with  forward  speed.  Furthermore,  the  cyclic  pitch  is  adjusted  to 
suppress  first  harmonic  cyclic  flapping  and,  therefore, 

to  eliminate  rotor  hub  moments. 

The  following  auto-pilot  equations,  taken  from  reference  (l) 
represent  the  strategy  whereby  the  controls  (Qc  t  C  •  6  )  are  adjusted 
in  order  to  reach  the  desired  thrust  (T=T.  )  and  moments  (?K=?E0= 0, 
RK=RMc=0) .  The  cross-coupling,  B,  accounts  for  the  coupxng  between 
roll  and  pitch  in  a  helicopter  rotor. 


Zc  ft  r  £  *  Ac  ( rt  -  r ) 

t,  e>  *■  ft  »  A,  [i/Vj  -  PHe )  B(RM  -  RH c)J  ^ 

t  f  r  £  -  A,  [(.RM-'RHe)  *  6  IpH-PK)]  (9) 

The  final  form  of  the  auto-pilot  equations  are  found  from 
evaluation  of  the  instantaneous  thrust  and  moments  (T,PK,RF.)  in 
terms  of  the  flapping  angle, £  . 

6  =  -4  -  AX/3  ikic  (1°) 

ft,  -  'A  -  -  cv^'j/s 

rA  (ll) 

1  ~ , 


fl  -  -  _  JilSA-JD.  I  .£yLr-Li££*  V-  t  4 

c,  r,  r  ^  /'i 


~  A.ii 


(12) 


2.3.  method  of  Solution 

For  the  whole  system,  we  combine  all  the  above  equations, 
and  write  them  in  matrix  form,  £» 
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(14) 


The  following  blade  parameters  and  flight  conditions  are 
used  in  equations  (13)  and  (14)  for  the  numerical  examples  to 
follow. 


/*  =0.3 

<r  =5.0 
e  =o.o 5 

Z  =o.i 


■■ 
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Cr/C&=Q. 0050/0. 314=0. 0159 

ct/r<.i=o.o 

Gm/{T'A=0.0 
Pi  =0.0 
=0.2 

<^  =  |  I  (j/^-  £>' )*  '  /*.!}*= 8. 330 12x10  1 
8, =0.0 

tJ  =3-5  .  P=1 .12  ,  «/>  neglected 

We  select  ft, ,  Q  ,  C  ,  $  ,  p  ,  £> ,  (ft  1  £*»  and  as  the  state  variables 
and  rewrite  equations  (13)  and  (14)  in  first-order  state  variable 
form. 

Y(1)=Y(2) 

Y ( 2)  =  —  [1 . 2544+( 0 . 3cosx-0 . 2)  ( 0 .62>0 .  l8?5sinx)J  Y(  1) 
-(0.625*-0.l875sinx)Y(2) 

+( 0 . 62>0 . 1875sinx)  ( 1+0 . 3sinx)  Y(  3) 

+0.025(0. 62 5+0 . 18?  5sinx) Y ( 4) 

-(0 .00520632+0 .00156l89sinx) 

Y(3)=Y(4) 

Y(4)  =  -(0.62>0.l875sinx)(0.25)Y(l) 

-12.25Y(3) 

-(0.625*-0.l875sinx)(0.25)Y(4) 

+11.25Y(5) 

+11.25sinx  *Y(9)  f  (15) 

+11.25cosx  *Y(?) 

Y(5)=Y(6) 

Y(6)=-(AC/Tt )*0.25088Y( l)-Y(6)/Tt +Te (0.0159) 

Y(7)=Y(8) 

Y(8)=-(Aj/T,  )(0.05088)(0.40?0cosx+sinx)Y(l) 

-Y(8)/T, 

Y(9)=Y(10) 

Y(10)  =  -(a,/t/  ) ( 0 . 05088) ( 0 . 40?04sinx-cosx) Y ( l)  , 

-y(10)/t, 

where  x=  y 

Y(l)=  9 
Y(2)=.d 

y(3)=  e 

y(4)=  e 
Y(5)=  6, ( 


-  &  - 


y(6)=  e. 

Y(?)=  £ 

Y(8)=  <o 

Y(9)=  £ 

Y(10)=  £ 

We  solve  the  first-order  differential  equation  set  (15)  by  the 
Runge-Kutta  program  provided  by  IBM  Scientific  Subroutine  Package  (SSP) . 
Figure  2  shows  a  typical  response  of  a  stable  system  with  all  zero 
initial  conditions. 
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31_0PTIMALITY _C RITERI A_AKD  _S  EARCH  _F0R_0PTIMUM .POINTS 

3.1.  Optimality  Criteria 

For  a  system  to  be  optimal,  some  kinds  of  criteria  must  be  used 
to  decide  on  the  utility  of  the  solution. 

In  our  case,  for  a  stable  system,  the  angles  £  ,£  >  $i  will 
ultimately  reach  the  final,  stable  position.  It  is  our  desire  that 
these  final  values  are  reached  in  as  short  a  time  as  possible.  Thus, 
we  choose  as  a  cost  function  the  time  required  (  i.e.  the  number  of 
rotor  revolutions  required  )  for  all  controls  to  be  withir.  ±0.5*  of 
their  final  values.  To  do  this  we  designate  TTC  ,TTt,  and  TTj  as  the 
respective  times  for  ^  ,  and  £  to  converge  to  within  ±.0.5'  of 
a  final  value;  and  we  designate  the  largest  of  these  three  as  Tmax, 
Figure  2.  Then,  we  use  the  minimum  of  Tmax  as  the  optimality  criterion. 

3.2.  Search  for  Optimum  Points 

As  we  have  seen  in  previous  sections,  it  would  be  quite  diffi¬ 
cult  to  obtain  explicit  expression  for  TT0 ,  TTt ,  TT^  ,  and  Tmax. 
Furthermore,  analytical  derivatives  of  Tmax  with  respect  to  A„,  A,  , 

Te,  and  T(  are  not  readily  available,  so  we  use  the  derivative-free 
method  in  our  search  program.  That  is,  we  evaluate  functions  only; 
and  no  derivatives  are  involved. 

The  controller  used  here  has  four  parameters  that  must  be  con¬ 
sidered:  A at  A,,  T0,  T,  .  The  procedure  is  as  follows, 

(1)  A  base  point  is  chosen  and  Tmax  is  evaluated. 

(2)  Local  searches  are  made  by  stepping  Ac  a  distance  0.2  to  each 
side  (  Ml  direction  in  Table  1  )  and  by  evaluating  Tmax  to  see  if 
a  lower  Tmax  is  obtained, 

(3)  If  there  is  no  Tmax  decrease,  we  do  the  same  to  K2  direction, 
and  then  to  M3,  M4,  ...,  until  M42  direction,  or  until  a  decrease 
is  found,  see  Table  1.  We  select  0.2  as  increment  for  A(  ,  0.3;'  for 
T„  ,  0.2,7  for  T,  . 

(4)  If  there  is  a  Tmax  decrease,  we  then  use  the  new  point  as  a  base 
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Increments  of  Parameters 


Direction 

Ac 

A. 

T-„*  r- 

Ml 

0.2 

0.0 

0.0 

0.0 

M2 

0.0 

0.2 

0.0 

0.0 

Change  one  parame- 

M3 

0.0 

0.0 

0.3 

0.0 

ter  at  a  time. 

M4 

0.0 

0.0 

0.0 

0.2 

M5 

0.2 

0.2 

0.0 

0.0 

M6 

0.2 

-0.2 

0.0 

0.0 

M? 

0.2 

0.0 

0.3 

0.0 

M8 

0.2 

0.0 

-0.3 

0.0 

M9 

0.2 

0.0 

0.0 

0.2 

M10 

0.2 

0.0 

0.0 

-0.2 

Change  two  parame- 

Mil 

0.0 

0.2 

0.3 

0.0 

ters  at  a  time. 

M12 

0.0 

0.2 

-0.3 

0.0 

M13 

0.0 

0.2 

0.0 

0.2 

M14 

0.0 

0.2 

0.0 

-0.2 

M15 

0.0 

0.0 

0.3 

0.2 

Ml6 

0.0 

0.0 

0.3 

-0.2 

Ml? 

0.2 

0.2 

0.3 

0.0 

Ml  8 

0.2 

0.2 

-0.3 

0.0 

M19 

0.2 

-0.2 

0.3 

0.0 

M20 

-0.2 

0.2 

0.3 

0.0 

M21 

0.2 

0.2 

0.0 

0.2 

M22 

0.2 

0.2 

0.0 

-0.2 

M23 

0.2 

-0.2 

0.0 

0.2 

Change  three  para- 

M24 

-0.2 

0.2 

0.0 

0.2 

meters  at  a  time. 

M25 

0.2 

0.0 

0.3 

0.2 

M26 

0.2 

0.0 

0.3 

-0.2 

M2? 

0.2 

0.0 

-0.3 

0.2 

M28 

-0.2 

0.0 

0.3 

0.2 

M29 

0.0 

0.2 

0.3 

0.2 

M30 

0.0 

0.2 

0.3 

-0.2 

M31 

0.0 

0.2 

-0.3 

0.2 

M32 

0.0 

-0.2 

0.3 

0.2 

M33 

0.2 

0.2 

0.3 

0.2 

M34 

0.2 

0.2 

0.3 

-0.2 

M35 

0.2 

0.2 

-0.3 

0.2 

M36 

0.2 

-0.2 

0.3 

0.2 

M37 

-0.2 

0.2 

0.3 

0.2 

Change  four  para- 

M38 

0.2 

0.2 

-0.3 

-0.2 

meters  at  a  time. 

M39 

0.2 

-0.2 

0.3 

-0.2 

M40 

-0.2 

0.2 

0.3 

-0.2 

M41 

0.2 

-0.2 

-0.3 

0.2 

M42 

-0.2 

0.2 

-0.3 

0.2 

Table  1  Searching  Directions 


point  and  go  back  to  (2).  If  there  is  no  decrease  in  Tmax,  then  we 
assume  that  a  local  minimum  point  is  found. 

A  flow  diagram  illustrating  the  above  procedure  is  given  in 
Figure  3* 


Figure  3  Flow  Diagram 


4.1.  Stepping  Distance  ; 

Although,  in  principle,  a  smaller  stepping  distance  implies  a  more 
accurate  convergence,  from  an  engineering  point  of  view,  it  is  often  more 
efficient  to  take  larger  steps  such  that  a  meaningful  change  in  cost 
function  can  be  found.  Thus,  in  this  case,  we  have  used  steps  0.2  for  A,, 
and  A(  ,  0.3  71  for  T0  ,  and  0.2  71  for  T/  . 

4.2.  Choice  of  Stable  Points  as  Starting  Points 

For  the  starting  points,  the  system  must  be  stable,  otherwise  Tmax  i 

will  be  equal  to  infinity  (in  our  program  '.ax  is  approximately  equal 
to  40  cycles),  and  the  search  prog:  -  .m  .  tall. 

4.3.  Local  Minimum  and  Global  Mlr.lr 

Since  local  minima  are  potent  '  possible,  several  different 
starting  points  are  considered.  T mould  .nese  converge  to  several  different 
local  minima,  the  global  minimum  c an  be  chosen  from  the  local  ones. 

4.4.  Optimality  Criteria  , 

The  minimum  Tmax  is  selected  as  the  optimality  criterion  in  our 
case.  As  we  can  see  in  Figure  4  and  5,  however,  the  difference  in  Tmax 
between  two  local  optima  Is  only  1.26  cycles.  Or.  the  other  hard,  or.e 
optimum  has  larger  oscillations  (in  the  steady-state)  than  the  other. 

Therefore,  it  might  be  good  in  future  studies  to  consider  these  oscilla¬ 
tions  in  the  selection  of  optimality  criteria.  ( 

4.5.  Optimal  Control  Setting 

With  blade  parameters  and  flight  conditions  as  shown  ir.  pages 
7  and  8,  and  with  the  minimum  Tmax  criterion,  the  optimal  control  setting 


■(i  is  as  follow, 


0S  Degrees  0C  Degrees  0O  Degrees 


6 


mm '■3~ 


L 


TT*  =  4.92  cycles 
Tmax  =  TT£  =  5-29  cycles 
6,  («°  )  =  6.5627  degree 
^(0=)  =  1.1347  degree 
£{(<*)  =-4. 5841  degree 

The  optimum  results  (Figure  4)  can  be  compared  with  another  local 
optimal  (Figure  5).  For  the  optimum  settings,  in  Figure  4,  the  low  cyclic 
time  constant  (0.6  ft)  causes  ±0.5J  oscillations  which  are  on  the  boundary 
of  the  accepted  level.  In  Figure  5»  a  larger  time  constant  is  chosen  (T, 

=  1.2  71.).  This  reduces  the  oscillations  to  +0.1°,  but  also  necessitates 
lower  gains  (A^,  A,  reduced  from  2. 6, 3. 6  to  2.0,  2.0).  The  lower  gains 
imply  that  the  over-all  convergence  of  the  mean  is  slowed,  making  the  mear 
the  critical  criterion  for  convergence.  It  should  be  noted  that,  for  an 
analysis  with  more  than  one  blade  (we  only  have  one),  the  oscillations 
would  greatly  decrease  since  the  rotor  would  filter  out  or.ce-per-rev 
from  the  controller.  This  would  alter  the  optimum  in  Figure  4  (were  ±0.5° 
is  critical)  but  would  not  greatly  charge  the  optimum  in  Figure  5  since 
the  oscillations  are  not  driving  the  solution. 

In  conclusion,  we  ca.n  say  that  the  gains  and  time  constants  found 
in  Reference  (l)  are  very  close  to  the  optimum  found  here.  On  the  other 
hand,  the  research  in  this  report  shows  that  there  are  two  local  optima 
with  nearly  equal  convergence.  A  comparison  of  these  indicates  that  ar. 
improved  optimality  criterion  might  be  obtained  from  a  more  stringent 
penalty  on  steady-state  oscillations. 
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6.  NOMENCLATURE 


c 

c 

C, 

cL 

Cm 

Sn 

Cs 

S 

CT 

S 

d 

e 

h 


area  of  the  blade  section 
slope  of  the  lift  curve 
gains 
coupling 
blade  chord,  m 

nondimensionali zed  blade  chord, 
blade  lift  coefficient 
normalized  blade  moment  coefficient, 
blade  moment  coefficient 
normalized  blade  moment  coefficient, 

root  moment  coefficient 


_i_ 

<r  a 


m 

t  a 


normalized  root  moment  coefficient 

blade  thrust  coefficient 
normalized  blade  thrust  ccefficien 


C 


_s 
<r  a 


t; , 


vTa 


length  of  blade  section 
inertial  ratio, 
lift  force 


Jr. 


1^ 


nondimens iona  1  i zed  lift  fur..e 


in 


Sl'lh 

constant  term  in  algorithmic  f<ua  Kut  t.t 


method 


total  inertia,  I  +  mr‘ 

y 


blade  section  inertias,  kg-m“ 

total  number  of  periods  of  integration 


>•  » 


ul  Cl 
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VK1 


reduced  frequency,  y  ^ 

rate  gains,  such  that  ^  is  number  of 
radians  to  full  application  of  linear 
control  predictions 

spring  stiffness  at  the  center  of  rotation 
flapping  spring  constant 
torsional  spring  constant 


number  of  blades 


lift  on  blade 

rolling  moment  at  hub 

mass  of  the  blade  section,  number  of 

controls 

pitching  moment  at  hub 
moment  about  the  pitch  axis 

nondimensionaliced  moment  about  the 

Mg 

pitch,  — 

number  of  second-order  degrees  of  freedom 

dimensionless  rotating  flapping  frequency 

slope  of  the  pi  ten  moment  coefficient 

pitch  moment  desired 

roll  moment  desired 

blade  radial  coordinate 

steady  root  moment  on  hub  in  rotating 

sv stem 


J 


t 


time 
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T 

U 


14'ViTz 

V 

IT 

ci 

dc 

oc. 

C(x'0(y  <c(z 
/* 

/^o'/^s  '/^c 

r 

9 

©e 

Qo'©c'©s 

A 

P 

or 

To'Ti 


thrust  on  the  blade 

total  velocity  of  blade  section  relative 
to  air 

vertical  component  of  air  speed 
horizontal  component  of  air  speed 
components  of  T/  at  the  blade's  o.g. 
vertical  component  of  wind  speed 
velocity  of  the  blade's  center  of  gravity 
angle  of  attack 
critical  angle  of  attack 
angular  velocity  of  blade  at  center 
of  gravity 
components  of  c i 

flapping  angle,  rad.  ,  j3Q  *  jSs  si  n  ^ 

+  pc  cos  y 

components  of  flapping 
lock  number, 

Iy+mr 

total  ‘pitch  angle,  rad,  0^  +  0  ©s  sin\h 
+  0c  cos  y- 

elastic  portion  of  pitch  angle 
collective  and  cyclic  pitch 
inflow  ratio 
advance  ratio, 

flr 

density  of  air 

rotor  solidity,  J?- 
/l  r 


- 


time  constants 
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<p 

inflow 

angle 

¥ 

rotor 

azimuth  angle. 

f=fit 

SI 

angular  speed  at  the 

axes  of 

(J& 

pi  tch 

frequency 

(•) 

!-  (  ) 
dt  1  1 

(*) 

—  (  ) 
d>  1 

=  ~  --  (  ) 
ndt  '  J 
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